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OVERVIEW 

• MOTIVATION - TURBINE ENDWALL HEAT TRANSFER 

• DESCRIPTION OF NEW MODEL 

• RESULTS 

1. FLAT PLATE 

2. ANNULAR TURBINE CASCADE 

3. TURBINE ENDWALL HEAT TRANSFER 

4. SUPERSONIC COMPRESSOR BLADE 

• SUMMARY 



EXPERIMENTAL ENDWALL STANTON NUMBER CONTOURS 
AS A FUNCTION OF 5 inlet AND Re chord 



RVC3D (ROTOR VTSCOUS C0P E_3.J>) 

BY R. V. CHIMA 
DESCRIPTION 

. EULER OR NAVIER-STOKES ANALYSIS 
FOR STEADY 3-D FLOWS IN TURBOMACHINERY 

FEATURES 

• CARTESIAN FORMULATION, ROTATION ABOUT X-AXIS 
'RECTANGULAR OR ANNULAR GEOMETRIES 

• SOLVES NAVIER-STOKES EQUATIONS 

THIN-LAYER FORMULATION, (NO STREAMWISE VISCOUS TERMS) 
RETAINS HUB-TO-TIP & BLADE-TO-BLADE VISCOUS TERMS 
BALDWIN-LOMAX OR CEBECI-SMITH TURBULENCE MODEL 
SIMPLE TIP CLEARANCE MODEL 

• NODE-CENTERED FINITE-DIFFERENCE FORMULATION 
EXPLICIT 4-STAGE RUNGE-KUTTA TIME-MARCHING SCHEME 
2ND + 4TH ORDER ARTIFICIAL VISCOSITY, EIGENVALUE SCALING 
VARIABLE AUj & IMPLICIT RESIDUAL SMOOTHING 

HIGHLY VECTORIZED & AUTOTASKED FOR CRAY Y-MP 

• STACKED C-TYPE GRIDS 







TURBULENT VISCOSITY PROFILE 



INTERMITTENT 

TURBULENCE 

OUTER LAYER 
INNER LAYER 


rF,RLCT-SMTTH &: BALDW TN-LOM AX MODELS 

TNNER LAVER: PRANDTL-VA N DRIEST FORMULATION 

CEBECI-SMITH 

Hi = pl 2 \du/dy\ 
l = nyD 

D = l- exp(—y + /A + ) VAN DRIEST DAMPING 

BALDWIN-LOMAX 

Hi = P^|w| 


OUTER. T.AVER.; CLAUSE * FORMULATION 


CEBECI-SMITH 


Ho 
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Kfr/6'u c 
i + S.5 (|) 
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KLEBANOFF INTERMITTEN CY FUNCTION 


BALDWIN-LOMAX 

ys r* J Vmoz/mar 

P° ~ PI w | wake option 

f(y) = y\“\ D 


< 



P A T,PWTN-LOMAX MODEL ANALYSIS 

(SEE PAPER FOR DETAILS) 

1. ASSUME SUBLAYER-WALL-WAKE VELOCITY PROFILE 

2. CALCULATE BALDWIN-LOMAX FUNCTION /(») 

MAX. OCCURS AT y m „ = .6465 
INDEPENDENT OF PRESURE GRADIENT 

NO MAX. FOR INFINITELY FAVORABLE dp/dx 

3. SPURIOUS MAX. CAN OCCUR AT EDGE OF VISCOUS SUBLAYER 
MOST LIKELY AT LOW Re & FAVORABLE dp/dx 



1. VELOCITY PROFILES, Re* = 105,000 




2. B-L FUNCTION /(„) 3. SPURIOUS MAXIMUM IN f(y) 


PROPOSED TURBULENCE MODEL 

TNNF-R LAYER (SIMILAR TO BALDWIN-LOMAX) 

Pi = pi 2 M 
l = tzyD 

D = 1 - exp(— y + /A + ) 
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PROPOSED TTIPRITT/ENOE MODEL 


PRESSURE GRADIENT EFFECTS 

• ACCELERATING FLOWS TEND TO REL AMIN ARIZE 

• MODELLED BY INCREASING A+ IN FAVORABLE dp/ds 

• CEBECI’S EXPRESSION FOR A + USED: 


A * = 


26 

y'l + 11.8p+ 


P + 


v dp 
pu* 3 ds 


• PRESSURE GRADIENT EVALUATED USING: 


0p_ 

ds ~ |V.| 


•Vp 


• “EDGE VELOCITY” V e EVALUATED AT A GRID LINE FAR ENOUGH 
FROM THE WALL TO GIVE THE GENERAL FLOW DIRECTION 

• KAYS-MOFFATT EXPRESSION WAS TESTED, EFFECTS TOO STRONG 


PROPOSED TURBULENCE MODEL 

LOCAL SHEAR MODEL 

• IN STRONGLY ACCELERATING FLOWS r + DECREASES WITH y + 

• MODELLED B"Y R.EPLA.CII^G Twait WITH r(y) IN D 

D = 1 - exp(-y + /A + ) 


• ERROR IN ORIGINAL PAPER - USED p, |w| ONLY 

• USED BY KAYS, PATANKAR-SPALDING, OTHERS 

• ALSO USED TO AVOID PROBLEMS AT SEPARATION WHEN r wa „ -* 0 


PROPOSED TURBULENCE MODEL 


OUTER LAYER 


Vo 


7 

Cwk 

Cxub 


Kp~f min | c wk y(\V maz \ - \V min \) 

0.825 - 

0.55 


PPOPOSED TURBULENC E MOPEL 

nTTTF.R T/ AVER. - FUNCTI QN_Z 
. DEFINE F = sfdy 

• INTEGRATE BY PARTS ASSUMING |w| -» 0 AS y — 6 


F = y\u\dy 



= «y|o - lo udy 

= - u)dy 

F = 6’u e 

. USE F DIRECTLY IN CEBECI-SMITH OUTER FORMULATION 

• ELIMINATES CONSTANT 

• DOES NOT REQUIRE KNOWLEDGE OF 6 C)R u c 

. DISCOVERED INDEPENDENTLY BY D. A. JOHNSON, AIAA 92-0020 


• y IS THE CENTROID OF THE /(y) CURVE 

J*f(y)dy = jl f(y) d y 

• EVALUATE USING COLE’S VELOCITY PROFILES 

n y/6 
0 .5 

.5 .55 
oo .606 

• USE EQUILIBRIUM VALUE Cxub = y/ 6 = 55 



• LOWER OPTION IS A CONVENTIONAL WAKE MODEL 

. EVALUATE C wt BY EQUATING TWO OPTIONS, ASSUMING 

y sep = .606 S 
F sep = u e 6/2 
AV/u e « 1 

• GIVES = 0.825 
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PROPOSED TURBULENCE MODEL 

3-D IMPLEMENTATION 

• GRANVILLE BLENDING FUNCTION 

Vtff = Vo tanh — 

Vo 

• MODEL APPLIED INDEPENDENTLY IN BLADE-TO-BLADE (r?) 
AND SPANWISE «) DIRECTIONS 

. INNER LAYER - USE BULEEV LENGTH SCALE 

_ 2s,,s i: 

y ' ~ s, + s f + + si 

• OUTER LAYER - USE ACTUAL DISTANCE ACROSS PROFILE 

y 0 = s n OR S( 

• BLEND r, AND C PROFILES VECTORALLY 

Vturb = \JVU] Vt( 


FLAT PLATE VELOCITY PROFILES 



COMPARISON OF FLAT PLATE VELOCITY PROFILES 
TO SPALDING’S COMPOSITE LAW OF THE WALL 




STATIC PRESSURE DISTRIBUTION 
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FRACTION OF AXIAL CHORD 


COMPUTED & MEASURED PRESSURE DISTRIBUTIONS 
FOR THE ANNULAR TURBINE CASCADE 




COMPUTED & MEASURED LOSS COEFFICIENT PROFILES 
FOR THE ANNULAR TURBINE CASCADE 



COMPUTED & MEASURED EFFICIENCY CONTOURS 
IN THE WAKE OF THE ANNULAR TURBINE CASCADE 


TURBINE CASCADE 
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ORIGINAL PAGE IS 
OF POOR QUALITY 








TURBULENT UISCOSITY 
NEU WODEL 


CONTOUR LEUELS 

t. nonon 

2.000017 

:i . unarm 
■•I . unarm 

S.OUDflU 

G. onmiu 
t . unarm 

u. unurm 

H. imiMiii 
iimiiuinh 
ru.mmim 
vn. mini in 
nn.nnnnn 
^n.nonnn 
su.tiuuun 

f*H. iiRURU 



OWQtNAL PAGE !S 
OF POOR QUALITY 


SUMMARY 


- MOST LIKELY AT LOW Re AND FAVORABLE dp/ds 

• NEW TURBULENCE MODEL PROPOSED 

- INTEGRAL RELATIONS FOR AND i USED WITH OS MODEL 

- EFFECTS OF dp/ds MODELED 

- WAKE MODEL PROPOSED 

• FLAT PLATE 

- R-L & NEW MODEL AGREE WITH LAW OF THE WALL 

_ LOCAL SHEAR MOD. DOES NOT AGREE WITH LAW OF THE WAL 

• annular turbine 

- GOOD AGREEMENT WITH EXPT. PRESSURE DISTRIBUTION 

- WAKE MIXING UNDER-PREDICTED 

• TURBINE ENDWALL HEAT TRANSFER 

- VARIATIONS IN ENDWALL St WITH Re PREDICTED WELL 

- EFFECTS OF dp/ds IMPORTANT 

• TRANSONIC FAN 

- SHEAR LAYER FROM BOW SHOCK ACTS LIKE VISCOUS LAYER 

- new MODEL OVERPREDICTS L.E. im 

- b-l model predicts REASONABLE L.E. Pt 


